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The  propa.^.ation  of  •■.^'-.a'-:   sriocks    In  a   stratified   atraosphere 

by 

G,   3,   'ihitham 

Introduction 

In  a  pi-cvious  report  [l],  a  coinpietc  theory  vjas  presented 
of  the  propagation  of  weak  spherically  synmotric  shocks  In  air 
of  liniform  density,  pressure  etc:   Later,  the  theory  was  extended 
[2]  to  a  case  in  which  the  undisturbed  fluid  into  vjhich  the  shock 
moved  vjas  not  uniform,  but  the  property  of  spherical  syinrnetry 
was  preserved.   In  the  present  report,  a  coi'ibination  of  the  meth- 
ods of  [l]  and  [2]  is  applied  to  the  problem  of  a  shock  (produced 
by  an  explosion  say)  propagating  in  the  atmosphere,  taking  into 
accoimt  the  variation  of  the  density,  pressure  etc.  of  the  air 
vjith  height  above  the  ground.   Since  the  air  is  stratified  and, 
in  the  most  important  problem,  the  shock  is  initially  spherical 
(approximately)  only  the  s^Trimctry  about  a  vertical  axis  remains; 
the  physical  quantities  are  novj  functions  of  two  space  variables 
and  the  time.   The  introduction  of  an  additional  ii-dependent  var- 
iable into  the  problem  Is  far  from  a  trivial  extension,  and  if 
the  theory  can  be  v:orkea  out  successfully  in  this  particular 
case,  4-t  vjill  st.rvc  as  a  prototype  for  many  other  investigations 
(in  the  theory  of  tliree  dimensional  steady  supersonic  flovj,  for 
example)  .   The  present  accoiont  constitutes  a  "progress  report" 
since  certain  points  require  further  consideration  and  justifica- 
tion before  a  final  report  can  be  written. 

Only  vjeak  disturbances  are  treated  here,  and  as  before  the 
linearized  (acoustic)  theory  forms  the  starting  point  of  the  method. 
The  linear  results  are  inadequate  in  certain  vjays,  as  explained  in 
[ij,  but  their  shortcomings  can  be  remedied  to  give  a  valid  descrip- 
tion of  the  flow.   The  v.-ork,  therefore,  falls  into  tv;o  sections:  (I) 
the  derivation  and  consideration  of  the  appropriate  results  of  linear 
theory,  and  (II)  the  Improvement  of  those  results  and  the  determina- 
tion of  the  shock.   The  first  section,  strangely  enough,  is  the  more 
troublesome  and  the  inore  inconplote ,  since  explicit  solutions  of  the 
linear  equations  are  not  Icnowr.  in  general,  and  we  have  to  be  satis- 
fied vrith  a  certain  expansion  solution.   This  expansion  , 
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hovjever,  is  of  considerable  interest  for  its  own  sake  and  is  of 
general  application  to  linear  hyperbolic  equations  of  second  order. 
The  first  section  extends  the  methods  of  treating  the  linear  theory, 
vjhich  vjere  developed  in  [2j,  to  include  an  '^.dditional  space  variable 
(the  treatment  of  linear  theory  is  trivial  in  the  problem  of  [l], 
since  the  explicit  solutiin  of  spherical  sound  x^faves  is  well-known), 
v/hile  the  second  section  extends  the  applicability  of  the  basic 
ideas  described  in  [ij. 

I.   The  linf.arizcd  theory 

The  physical  quantities  are  functions  of  the  distance,  r,  from 
the  (vertical)  axis  of  syi^etry,  the  height,  z,  above  the  ground, 
and  the  time  t.   If  u  and  w  are  the  velocity  coiaponents  in  the  di- 
rections of  the  r  and  z  ax^.s,  respectively,  the  density  is  p  and 
the  prv;ssure  is  p,  then  the  equations  of  momentujn  (including  the  ac- 
celeration of  gravity  g)  are 

(1)  u,  +  uu   +  wu  =  —  p 

^  t     r     z     p  -  r 

(2)  w.  +  uvj   +  WW  =  -—  p„-g, 
^    '  t  r     z    p  -^z  •=" 

and  the  equation  of  continuity  is 

(3)  Pt  +  up  I.  +  wp^  +  P(u^  +  p  +  w^)  =  0. 

It  may  be  assumed,  sincj  entrap;/  changes  are  negligible  for  weak 
disturbances,  that  p  is  a  given  function  of  p.   Then,  introducing 
the  velocity  potential  '-]),  equations  (l)  and  (2)  may  be  Integrated 
to  give  Bernoulli's  eouatlon 

(I4-)        i^  +  2^^r   "*"  ^z^  "t/  P  "*"  SZ  =  constant. 

In   equilibriuri  u  =  vj   =   0,    and   p   =    P(z),    p   =  R(z),    a  =  /dp/ dp  = 
A(z),    say, 

and   from    (2)    this    equilibriiyri  conrigiir-ation   is    determined   by 

(5)  P'(z)    +  gR(z)    =  0. 

To  obtain  the  linearized  equation  of  motion,  tho  deviations  of  <}),p, 
p,  am.d  a  from,  their  equilibrium^  values  (the  undisturbed  value 
of  4"  is  zc.ro)  are  asstuned  to  be  small,  and  in  the  above  equations 
only  the  first  order  terms  in  these  small  quantities  are  retained. 
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Equation    (I4.)    becoraes 

(6)  p  -  P(z)  =  -R(z)  ([)  t, 

and  (3)  becomes 

(7)      ^  +  i  +  ii:  +  i  ^(Ri  )  =  0. 

In  addition,  since  a   =  dp/dp,  p^  ^=^(p-P)^A  'z);  hence,  fro 
(6),  (7)  is 

(8)         ^tt  =  '^^-^     ['^rr  -^  Ftr  ^  ^zz  "^  lUf  ''^^ 

The  solution  of  this  equation  for  'f*  is  not  kn^^wn,  but  an  ex- 
pansion valid  near  the  fr.int  .f  the  disturbance  will  prove  to 
be  of  value.   To  motivate  this  expansion  we  consider  tvjo  special 
cases:  (i)  the  case  of  a  uniform  atinosphcre  with  A  equal  to  a  con- 
stant Ji     and  -^   a  function  :•£   r  alone  (this  is  the  problem  of  cy» 
lindrical  sound  vjaves)  ,  and  (ii)  the  case  of  plane  waves  propagat- 
ing vertically  so  that   4"   i^  ^   function  only  of  z.   In  each  of 
these  problems  only  one  space  variable  appears,  and  the  equation 
reduces  to  the  type  considered  in  [2];  (  i)  and  (il)  are  no'":   con- 
sidered in  turn. 

( i)  Cylindrical  sound  waves 

In  this  case  (8)  reduces  to 

(9)  ^     +  ik  =  itt 

^^'  ^rr    r    ^^» 

o 

and  it  is  vrell-kn  ?v;n  that  the  solution  which  re^jrcsents  an  out- 
going; wave  starting  at  t  =  0  is       /. 

•'* ■  ^/" o       h(t' )  d t » 


(10)        4'  =  K  (pr/A^)  h(t)  =  J        /-     2 


I  V(t.tr)^  -  v^/l? 

Whoro  h(t)  is  an  arbitrary  function,  K   is  the  Besscl  function,  and 

o  '   ' 

p  denotes  the  operator  in  the  Hcaviside  calculus.   (Tho  pressure 
has  also  been  den  )ted  by  p,  but  no   confusion  will  ensue,  since  the 
pressure  does  not  appear  a.vcain  in  this  section)  . 

It  is  found  in  these  pr  ;ble:ns  (see  [l]  and  [2])  that  in  order 
to  determine  the  head  shock  the  linearized  solution  is  required  in 
the  region  v/here  A^/r  is  small,  here,  E,    is  the  characteristic  vari- 
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able t-r/A   which  mc-isurcs  the  time  after  the  nrriviil  of  the  first 
disturbance  (which  travels  on  ^  =  0)  .   Although  small  E,   corresponds 
to  the  front  of  the  Usturbance  in  the  linearize;'',  theory,  the  solu- 
tion for  sraall  E,   is  n'-)t  sufficient.   For  in  the  imor  rv'c  :1  (non-lin- 
ear) theory,  characteristics  are  continually  meeting  the  shock  so 
that  eventually  at  large  distances  the  appropriate  value  of  the 
characteristic  variable  near  the  shock  is  no  longer  sraall.   Hence, 
roughly  speaking,  we  require  the  solution  for  small  ^  and,  in  aildi- 
tion,  the  solution  for  large  r|  more  precisely,  the  solutions  for 
small  A4/^  is  needed. 

The  expansion  can  be  obtained  from  the  integral  In  (10),  hut 
is  found  more  quickly  from  the  Heaviside  representation  of  the  so- 
lution.  Small  values  of  E,   correspond  to  large  values  of  p  in  the 
Heaviside  representationj  in  fact,  l/p  plays  a  role  corresponding 
to  4  so  that  the  expansion  of  ^    for  small  values  of  k^E,/v   is  found 
fr.^im  the  Heaviside  representation  by  assuming  k,/-pv   is  small.   Ex- 
panding the  Bessel  function  for  large  values  of  its  argument, 
pr/A,  ,  and  interpreting  term  by  term,  we  have 


o- 


^      --P^A.     ].         Ill  ^9 


.2 


i\ 


I'  {  ^  p   r  ) 

=    (^)2   f(t-r/A^)    -  |(^)2   f,(t-r/A^)    +  ^Cf)^  Ut-:^  -  , . .  , 

where  

f(t)    =  JfrUt)    =i      \      Ht')it' 

and  f  (t)  denotes  the  nth  rcpoatcd  integral  from  0  to  t.   To  a 
first  order  approximation  only  the  first  term,  in  (11)  need  be  re- 
tained; using  it,  the  linearized  theory  may  be  Improved  and  the 
head  shock  determined  by  the  procedure  described  in  [ij.   The  de- 
tails are  not  given  here,  since  this  special  case  is  -nly  used  to 
m.^tivate  the  later  treatm.ent  :f  the  linear  equation  (8);  in  any 
case,  the  results  are  very  similar  to  these  of  [3]. 

( ii)  Vortical  propagation  of  plane  w.aves 

It  is  assumed  n )W  that  the  air  is  polytropic  i.e.  the  pressure 
is  equal  to  kp*^,  where  k  and  y  '-^'^    constants.   It  follows  then  from 
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(5)    uhat  R{z)    cc(z   -z)y-T,    and   t  lie  mi  fore   that 


(12) 


A  =  A^(z^-z)?, 


where  z  and  A_  are  constant';,   ■'>-j 
o      1 


z  ,  R  and  A  tend  to  zero,  but 
o' 


of  course  the  simple  polytrcvo.c  ceases  to  apply  to  the  atmosphere  be- 

2 

fere  z  =  z   Is  reached  (for  cjcamvle,  the  tropopause,  at  which  T  oc  A 

Is  discontinuous,  intervenes)  j  hc-jover,  the  behaviour  according  to  (12) 
of  A  near  z  =  z   will  prove  interesting  mathematical  Ijr,   VJith  these 
valties  of  R  and  A,  and  J)  a  fionction  of  z  alone,  (8)  rediices  to 


(13) 


''^tt  =  ^l^^o-^^^ 


I 


zz 


1    (i;Z 

Y-1  z^r^ 


It  is  convenient  to  introduce 

(11^) 

in  place  of  zj  then 


(IS) 


^  1 


where  v  =  p-C  3-y)/(  Y'l)  j  ^-^^  the  characteristics  are  t  +_/p   =  constant. 

I  2 1 

Replacing  (P^.  by  the  operational  form  p  d).  the  general  solution  of  (15) 


la 


(16) 


A  =Y'^    ['v4  (^P)l-^l(-^)  +  K,^4(f  )h2(t)j. 


whore  I  ,1  and  K^A  denote  Bessol  functions  of  Imaginary  argument,  and 

V  I  '<5'        V  I  H 

h^(t),hp(t)  arbitrarjr  function:j.   For  simplicity,  consider  the  special 
ca?:e  v  =  integer  f -r  which 

Then,  it  is  easy  to  pick  out  the  solution  representing  propagation  on 
the  'out£;olng'  characteristics  t+J>j7=  constant  (note:  with  the  definition 
(lli-),^  decreases  as  z  increases,  hence  the  outgoing  wave  hr- s  the  posi- 
tive sign  in  the  characterls bie  equation);  it  is 

(17)      cL  =:^''L(^p)j(t), 

d  ^ 

Interpreting  (1?),  we  have  the  solution 
(16)       I   =- 


Y 
where  L(^)  =  Y^      {rrr^)^  [r^-) ,    and  o(t)  is  an  arbitrai'y  function. 


J,  =Y^j(t.^)  .  viv«)^v-l.^(^,^,,__^,i^i|^l^.^( 
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where  j  (t)  denotes  the  nth  integral  of  j(t),  the  lower  limits 
being  the  value  of  t  +  ^"  at  the  head  of  the  wave.   Near  the 
head  of  the  wave,  the  successive  terms  decrease  in  magnitude 
provided  ^  4  Q»  ^^^  the  solution  is  of  exactly  the  same  general 
form  as  (ll)j  that  is, 

(19)  ^^^^^''hiO  +  ^'^^h^{0  +  ^.*, 

whore  ^   is  the  characteristic  variable  defined  to  be  zero  at  the 
head  of  the  wave,  f(^)  is  on  arbitr'iry  function  (j(t+/^)  =  f(5)), 
and  the  (j) '  s  are  certain  functions  of  the  space  variable.   In  the 
particular  case  (l8),  there  arc  only  a  finite  number  of  terras, 
but  for  non-integral  values  of  v  the  series  would  be  infinite. 
We  note,  however,  that  in  this  case  the  scries  is  valid  only 
for  small  ^  and  V  ^   0;    as^'-s»  0,  it  is  not  uniforroly  valid. 
For  fixed  g,  the  last  term  of  ( l6 )  dominates  for  sufficiently 
small  values  of  ?7   j  the  order  of  the  terns  should  be  reversed, 
and  the  expression  (l3)  {which  in  general  is  replaced  by  an  in- 
finite scries)  would  invjlve  successive  derivatives  of  j  (t). 
Thus,  vie  see  that  although  (19)  is  valid  in  the  initial  stages 
of  the  propagation,  it  miy  not  apply  in  l-'.tor  stages. 

( i i i )  Discussion  and  genoralization  jf   the  results  cf  (i)  and  (ii) 

An  expansion  ''f  the  form  (19)  can  be  found  quite  generally 
for  an  equation  of  the  form 

(20)  (j)^^  =  /,2(z)  ji^^  +  B(z)  f.^  +  C(z)^  j  , 

even  though  the  explicit  solution  is  not  known,  by  substituting 
(19)  in  the  equation  and  equating  the  coefficients  of  the  f ' s 
to  zero.   This  was  used  extensively  in  [2j;  in  all  cas-js  (apart 
from  the  occurrence  of  singularities  as  in  the  problem  discussed 
above),  (19)  is  valid  for  small  5,  and  in  certain  cases  such  as 
(11)  its  v^.lidity  extends  into  the  region  ^^^/r   small.   The  ex- 
pansion is  intimately  connected  v;ith  the  W.K.B.  method  :f   find- 
ing asymptotic  solutions  of  ordinary  dif fer._ntial  equations. 
For,  the  operational  form  -f  (20)  is 

'I'zz  "^  ^(^^•i'z  ""f  (^^  -  P^A^(2)j.|.  =  0, 

and   the   asyiTiptotic    solution  for   large   p    (corresponding   to    sm_all 
4)    is   obtained   in   the   V/.K.B.   method   bj  the    substitution 
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*-e-P^(^)   (  ]S-hz)    ^  ^l-^^hz)    ^    ...!  f(p), 

whcro   f(p)   is  nn  arbitrnry  function  of  p.   It  is  founri  that 
q(z)  =Ja"  dz,  hence,  letting  f(p)  be  the  representation  of 
f(t),  the  Gxpansion  is  exactly  (19). 

Another  interesting  fuatur^  of  ( IQ )  is  in  the  significance 
of  '|)^"  .   Suppose,  there  is  a  :Uscontinuity  in  f(^),  then,  since 
f^(.^)  is  continu  ;us  for  n  >  1,  (19)  shows  that  an  initial  dis- 
continuity "would  be  propagated  along  a  characteristic  and  its 
magnitude  is  proportional  to  '|j  ^ ''  .   (Since  <\)    is  the  velocity 
potential,  discontinuities  of  >-[;  itself  W':uld  not  arise.   How- 
ever, (j)i^  is  proportional  to  the  pressure  Increase,  and 
'|;^  f-o  (|)(°^ft(^)  +  (},(^)f(^)  +  ...  .  therefore,  the,  same  result 
is  obtained  fir  pressure  jumps  by  considering  a  discontinuity 
in  f'(^)).   On  the  other  hand,  the  propagation  of  -  such  discon- 
tinuities can  bo  deduced  by  certain  standard  methods  directly 
from  the  coefficients  oT   oquati;n  (20);  if  this  Is  done,  exactly 
the  same  function  f{)^'^'  is  found  f.)r  the  magnitude  of  the  dis- 
continuity as  it  moves  *long  the  characteristic.   Thus  (19)  al- 
ready contains  the  correct  (linearized)  propagation  of  discon- 
tinuities.  Moreover,  since  we  have  a  complete  expansion  near  the 
head  of  the  wave,  additional  Inf  .rmation  is  furnished  about  the 
variation  in  the  magnitude  of  the  higher  derlvativos.   For  v.x- 
ample,  suppose  as  a  special  case  ->f  (19)  wo  have 

Then,  the  discontinuities  in  <\-,     and  h      on  ^  =  0  are  proportional 
to  ^\>  and  -  A"  4   ',  'the  discontinuities  in  ^+-+-5  ^^^^    <-ii^'  i^z  °^ 

^  =  0  are  proportional  to  cj)^^'',  dcb^^'^/dz  -  A'-^cj)^^^,  and 

.|A~^!J)^^^  -  2A"-^dc|)^°V^'z  +  A"^.l)^"^^-U/dzV  ,  respectively  and  so  on. 

^  In.lependcntly,  M.  Kline  has  c  nsidered  a  special  case  of 

(19)  (and  of  the  extensions  \Jhich  will  be  described,  in  (iv))j 
he  conoi.^ers  periodic  waves  an;*,  thus   f(^)  =  o"^  ^.      Although 
it  is  convenient  in  our  vrork  to  take  the  lovjer  limits  of  the 
integrals  in  (19)  as  ^  =  0,  the  expansion  Is  a  solution  of  the 
equation  vjhatever  the  limits  are,  or  even  if  the  Integrals  .ore 
indefinite;  hence,   f^(^)  ^^J   '^^  taken  as  e-^^^^/w^.   In  this 
way  Kline  has  an  expansion  valid  for  large  w  (for  any  E,)  , 
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( iv )  l!!b&   f;cneral  problem 

The  extension  of  ( 19 )  to  include  an  extra  space  variable  is 
now  discussed  and  in  a  similc.r  way,  any  nuraber  of  space  variables 
could  be  treated.   V/hen  an  additional  space  variable  is  introduced, 
bheie  no  longer  exists  a  unicp  e  set  of  characteristic  surfaces 
4(t,r,z)  =  constant,  since  the  characteristic  conoids  froin  any 
curve  in  the  plane  t  =  0  have  an  envelope  vjhich  is  a  characteris- 
tic surface.   Hoviever,  the  vjave  front  t  =  a(r,z)  is  uniquely  de- 
fined, ir\d   the  obvious  extension  of  the  4  used  in  (19)  is  $  ^ 
^  t  -  c.{r,z)  ,    so  that  at  any  point  in  the   (r,z)  plane,  £,   meas- 
ures the  tine  which  has  elapsed  since  the  passai^e  of  the  wave 
front.   In  addition,  the  i|)'s  of  (19)  now  becorae  fu.nc  lions  of  r  and 
z.   ^^s  before,  this  expansion  will  be  valid  near  the  head  of  the 
wave  since  f  /f   -,  =  0(§)  ,  and  in  corto.in  cases  the  i-'0';^ion  of  va- 
lidity mo-'j   be  extended  to  inclLide  points  at  lar!;ie  distances,  as 
in  (i),   ''^hen  the  polytropic  assujiiption  is  made  a  behaviour  simi- 
lar to  (ii)  is  expected;  in  fact,  near  A=0,  any  "wave  will  tend, 
by  refraction,  to  become  parall^.l  to  the  line  A=0.  hence  the  simi- 
larity to  (ii)  should  be  ^'^ery  close. 

By  substitution  of  the  extension  of  ( 19 )  into  the  equation 
(8) ,  we  have 

(21)  4^4   =   l/-^('^),? 

(22)  2a   A^'-^)  +  2a   .|.^^)-i-ia   +^+a   +  ^i^a   {^^'^  =  0 , 
^      '  r  ^r       z  ^z      ■  rr    r     zz    R(z)  zP^       ' 

(23)  2.,  *<."'  .  2a^  e^    *(y,   ^f   *  .,,  .  ^^^flhz}*'"' 

The  first  eqiiation  is  the  well-knoi>;n  eiconal  equation  vjhich 
must  be  satisfied  by  the  v;ave  front.   If,  for  cx<ari.pl.^ ,  the  initial 
vave  front  at  t  =  0  is  known  i.e.  the  solution  of  a(r,z)  =  0  is 
given,  (21)  is  sufficic^nt  to  determine  the  function  ^■-(r,z).   Geo- 
metrically, tho  surface   t  =  a(r,z)  in  the  (t,r,z)  space  is  the 
envelop.:,  of  the  characteristic  conoids  drax^jn  from  points  on  the 
initial  wave  front,  the  successive  wave  fronts  in  the  (r,z)  plane 
for  different  times  are  the  projections  (on  thu  (r,z)  plane)  of 
the  curves  of  intersection  of  the  suj?face  t  =  a(r,z)  with  planes 
t  constant. 
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It  is  convenient,  in  addition  to  the  set  of  viave  fronts  a  = 
constant,  to  consider  the  'rays'  which  are  the  orthogonal  tra- 
jectories of  the  vjave  fronts  (J^igure  1),   The  direction  cosines 
of  a  ray  at  any  point  are  (Aa  ,  A.a  )  (by  (21)  the  si^ii  of  the  squares 
is  1),  and  therefore  the  rate  of  change  of  any  quantity  along  a 
ray  is  determined  by  the  operator 


VJave  /fron 


(2k) 


/■  +  Aa  ^ 
2.  2 


In  particular,  da/ds  =  A(g,  +a,„)  =  l/A,  therefore  the  ivavo  front  ex- 
pands  out  along  the  rays  with  velocity  A.   In  terras  of  the  ray  deri- 
vative the  left  hand  sides  of  (22)  and  (23)  become  d(|)^^V'^^-S  +  K(j;^^\ 
where  X  is  a  known  function  of  r  an  1  z,  assuninp;  that  o.{v,z)    has 
already  boon  dcternined.   The  equs.tlons  for  the  (f) '  s  thus  become 
simply  first  order  ordinary  differential  oqufitlons  along  the  rays. 
Again,  the  propagation  of  discontinuities  in  the  solution  of 
equation  (8)  can  be  studied  directly  by  an  alternative  n;.thcd  and 
exactly  equations  (21)  and  (22)  arc  obtained  for  the  wave  fronts 
and  the  magnitude  of  the  discontiuity.   For  this  reason  the  treat- 
ment of  those  equations  -  the  theory  of  geometrical  optics  and 
acoustics  (see  [i|])  -  is  already  knov/n;  the  rcl^:vant  part  is  repoatec 
here.   Equation  (22)  can  be  written 


^r^  r^ 


(o) 


^R)  +  5t(^J 


(o)' 


rR) 


0; 


^z'  z^    ■  '     » 

integrating  this  expression  over  the  region  ABCD  (Fig.  1)  contained 


\ 


t  ■  ■■■ 

t 

i  / 


..^ 


■  / 


■ic;    ■-■■ 
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bet\JOGn   ajiy    tvjo   wave    fr.'nts    aril    tvjo   adjacent   rays,    tht.    'livorgence 
thoorem  gives 

(25)  /(Ia^+ma^),|,^  ■>    rnlo-^   0, 

where    the    integration   is   over   t]ie    boundary  of  ARCD,    an'    ,-^  ,ra   are 
the    lircct    :.n  cosines    of   the     .atv^^rl  n.-rn'-il.      Since    (Aa^,    ^a    )  are 
tiio    -i.irecti.-n   cosines    of   the    r-.yo,    -C  a   +rna      =   0    >n  /.C    and   BD, 

while  on  -.B  and  CD,  ^^a      +ma   =  /.(u,^  +  ^?  =  lA«   Th.,rof  irt^ ,  (25) 
gives  approximately 

[■\=^''-^   rxRAoA]i-[|'^  ■'  rR./^icrA]^  =  0, 

wherc/^>cr  -"'.tn -ttes  tb-e  length  of  the  small  segment  of  wove  fnnt 
cut  off  by  the  'rays,  and  subscripts  o  and  1  indicate  that  the 
quantity  j.s  evaluated  at  i-B   anl  CD  respectively.   Defining;  the 
expansion  rotio  E  =  lim  ^x/'2\-^,  we  have  j;'''^   rRE/A  is  con- 

o 
stant  on  each  ray,  and  in  particular  equals  its  value  at  the  ini- 
tial vjave  front.   In  this  vjay,  wo  may  consider   a{r,z)  and 

(1)^'   (r,z)  to  be  determine:!. 

( n^ 
Equations  (23)  for  the  later  functions  '[>  '  ,    ov^n  th';uo;h  they 

are  reduced  by  (2l).)  to  equations  -ilong  the  rays  are  complicated  to 

deal  with,   ^^n  investigation  -f  them,  w-iuld  be  of  interest  in  yiel'l- 

ino  additional  information  ab.^iut  the  region  of  validity  of  (19); 

also  these  later  terms  show,  in  r.   similar  way  to  that  described  in 

CL1.1 ),  h'iw  discontinuities  in  higher  derivatives  vary,  acc-^rding  to 

the  linear  theory. 

II   Improvements  of  the  linearize^,  theory  and  determinatiion  of 
the  head  shock. 

Aoc.-.rding  to  the  linearizt^d  thc^.ry,  in  the  early  stages  of 
the  propagation  and  even  later  in  certain  cases,  we  have  a  solu- 
tion )f  the  form  (19),  and  as  o  first  approxi^n.-.ti  _n  we  take  only 
the  first  term 

(26)  ,])  =  X(r,z)f(t-a(r,z)), 

where  fir  convenience  in  vjritm.,  \^''^'    has  been  replaced  by  X.   It 
is  n -w  sh-:wn  how  (26)  moy  be  im^yroved  and  the  head  sh)ck  deter- 
mined.  The  method  will  apply  wh..rever  |  is  of  the  apor-^ximate  form 
shown  in  (26)  even  if  it  no  longer  comes  from  an  exp.anslon  (19),  For 


-11- 

exnraple,  in  the  prjblem  describe--''  in  (ii),  (19)  did  not  apply  near 
V    -   0,    -ut,  near  ■*9  =  0,  another  expansion  vjh.ise  .'dominant  term 
was  still  of  the  f  :)rm  (26)  (with  different  X  and  f )  ,  wns  valid, 
therefore,  the  method  described  belov;  could  still  be  applie  1  in 
that  region. 

The  functirms  u.(r,z)  and  X(r,z)  are  assumed  to  be  known  from 
the  geometrical  acoustics  described  in  the  previous  section  and  f 
is  an  arbitrary  function  which  is  fixed  by  applying  some  appropri- 
ate boundary  or  initial  conditi -.n.   Unfortunately,  appropririte 
boundary  conditions  would  be  applied  in  regions  outside  the  range 
of  validity  of  the  expansion  (19).   For  example,  if  the  acoustic 
waves  in  I(i)  are  produced  by  the  small  motion  of  a  "cylindrical 
piston"  of  initial  radius  equal  to  zero,  the  boundary  con^lition  is 

applied  vrhere  r  is  small,  and  ;[/ ^'^^A  r~  I  h(t')dt'  there).   Hovjever, 
in  explosion  problems,  the  gen>-:.ral  form  of  f  is  known  from  examina- 
tion of  special  cases  such  as  the  problem  in  [lj«  therof 're,  a 
suitable  form  can  be  assumed  when  required. 

Henceforth,  we  shall  be  m   re  interested  in  tho  v.alues  of  u,w, 
and  a  thaji  in  the  velocity  potential  cf>j  these  are  n  iw  set  down. 
The  dominant  terms  in  u  and  vj  are 

(27)  u  =  ci^XF(t-a)  , 

(28)  w  =  a^XF(t-a)  , 

where  F(^)  =  -f ' (4) .   The  sound  speed  a  is  obtained  from  Bernoul- 
li's equation  (l|.),  but  it  involves  the  pressure-density  r(.lation; 
for  simplicity,  we  assuiae  thit  tiie  air  is  polytropic,  although  the 
general  cas:.  could  be  carrie  1  out  in  the  same  vjay.   Then 

(29)  a  ^  A  -  :gi  |,^  =  A  +  1:^  XF(t-a)  . 

In  ace  .rdance  with  the  hyp.,'bhesis  which  has  been  developed  and 
discussed  in  detail  in  previ  jus  ...apors  ([l],  [2],  [3]),  the  linear- 
ized theory,  which  becomes  ina-leqtiate  after  the  disturbance  has 
propagated  for  some  time,  may  be  corrected  quite  simply  to  provides 


For  integral  values  of  v,  tji.;  two  'expansions'  wer!_^  the  same 
expression  written  in  different  order;  from  the  present  point  of 
view,  hovrcver,  they  may  be  tr^^atod  as  distinct,  since  they  would  be 
in  general . 
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valid  description  of  the  motion.   In  the  exact  non-linear  theory, 
SiT-all  'wavelets'  (v/hOEe  paths  are  characteristic  surfaces  in  the 
(tjr^z)  space)  raov^.  at  the  local  sound  speed  a  relative  to  the  fluid 
VL^locity  (u,w).   The  linear  theory,  hov.'cver,  approxinatos  this  value 
by  the  undisturbed  sound  speed.  A,  on  the  grounds  th«it  a-A,  u,  w 
arc  small  compared  to  a.   But ,  as  a  wavelet  travels  the  error  intro- 
duced by  this  assumption  becoraes  large,  since  the  small  leviations 
in  tiie  .yroprigation  velocity  accuraulato .   To  correct  thie,  the  linear- 
iz-~.  characteristic  variable   t-a(r,z)  must  be  replaced  in  (26) -(29) 
^y  "C  (t,r,z),  whi  re  TT  is  a  sufficiently  accvrate  approxiraation  to  the 
exact  characteristic  variable.   Thus  (27),  (25),  (29)  become 

(30)  u  =  a^XFrcr), 

(31)  w  =  c^T^'ijZ^  , 

(32)  a  =  A  +  ^  y^{Tf), 

vjhcre  ■^(t,r,z)  remains  to  be  ■--.t.^rmined . 

A  surf  acc^  t,r,z)  =  constant  in  the  (t,r,z)  space  corrcspe^nds 
to  thv,  propagation  of  a  wavol^it,  and  gives  a  set  of  curves  in  the 
(r,z)  ph-ine ,  each  curve  corresponding  to  a  single  value  of  t  and  rep- 
rcSGntlTi.3  the  position  of  the  ijavulut  at  that  time.   It  is  easy  to 
show  that  the  velocity  of  the  v/avelet  normal  to  its^^lf  in  the  (r,z) 

plane  is  a  -  ZT  /  Vr.'2+T^2  .   xhe  vel'city  of  the  oarticlos  of  fluid 

r  z 
nor;nal  to  the  wavelet  is  (uTT  +  \nir) /  V^+T^    J  tnerefore  the  char- 
acteristic condition,  that  thi.-  v..l,'city  of  the  vjavolet  differs  from 


ty  b;\  the  soun: 

\    speed,  gi 

+  wrr 

^  1  n 

~ 

-  +  a. 

r    z   .    r 

z 

or 


(^^)  (--  -.u^^  +  w^-)2  =  a^fc^+TT^). 

Since  expressions  for  u,w  and  a  in  terms  of  x:  are  known,  (33)  ^-ives 
an  equatio:-]  to  determine  "c" .   Njvj,  al  though  TTdiffers  from  t-a  by  -n 
amount  which  becomes  large  (this  being  one  of  the  crucial  points  in 
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the    •-liscussi  m)  ,    the    :Ufforenci;    is    still    small   coraparud    to    a. 
Hence,   xire   may   Sc;t 

IZ   =   t    -    ^(r,z)    +   ix(i',z)  , 

where  [i/^    is  small.   Th' n,  subs  ti  tutinr-  In  (33)  ■'^nl  neglecting  terms 

of  second  order  in  small  quantities,  we  obtain 

ua   +  v:a       , 

(xising    (21)).      From    (30),    (31)    and    (32),    the   right  hand   side    of    (3lj.) 
is   ^( Y+l)XFtC)/A^.      Hence,    [i  may  be    taken   equal   to    p(r,z)P(^30,    vrhere 

(35)  ^   (3      +   a    S      =  -ilti  -4-    , 

from  (2l(.)  ,  the  left  hand  side  may  be  replaced  in  terms  -jf  a  ."'.eriva- 
tivo  along  a  ray,  to  give 

^^°'  ds    2   3 

To  define  p  uniquely,  it  is  stipulated  that  6=0  :n  the  initial 
wave  front  I  then  P  is  determined  by  integration  along  the  rays.   The 
solution  ( 30)  ,  (  31)  ,  ( 32)  is  njvj   completed  by  the  addition  of  the  Im- 
plicit  relation  f  or "^r-C  t ,  r,  z  )  : 

(37)  t  =  u(r,z)-|(Y+l)!3(r,z)F(c;)  +  T7. 

In  the  problems  of  principal  interest  (in  which  the  disturbance 
is  caused  by  an  explosion)  the  pressure  and  the  sound  speod  increase 
at  the  head  of  the  dlsturbancoj  thus,  PCO  is  initially  positive. 
It  foll.ws  from  (37)  that  th..  characteristic  surfaces  on  vJhich  F(^ 
Is  positive  vjill  eventually  intersect  the  characteristic  surfaces, 
t  -  u,  =  constant,  of  the  iindisturbed  region  ahead  where  FCW  =  0. 
(ivnother  way  -jf  expressing  the  same  thing  is  that,  in  the  (r,?) 
piano,  I'/avelets  carrying  greater  values  of  PfcJ  travel  faster  and 
overtake  the  others)  .   If  this  is  all-.iwe'l  to  occur,  the  solution  is 
no   longer  single  valued  an'  ceases  to  have  physical  significance; 
the  hua-";  shock  must  be  fitt-..  '.  in  this  region  in  ^rdor  that  the  char- 
acteristic surfaces  meet  it  and  are  cut  off  by  it  before  they  inter- 
sect each  other. 

i4.t  any  time,  the  fellowing  conditions  must  be  satisfied  across 


the  shock  in  the  (r,z)  pianos 

(a)  The  component  v,    of  the  particle  velocity  tan^^ential  to  the 
shock  must  be  continuous. 

(1:)  V  =  2(a-A)/(Y-l)  ivhoro  v   is  the  component  of  th^  particle  ve- 
locity normal  to  the  shock,  v   and  a  being  measure',  just  behin"  the 
shock. 
(c)  U  =  ^(A+a+v  )  vrhere  U   is  the  normal  velocity  of  the  shock. 

Txiose  are  deduced  from  the  Rankine-Hug'mist  shr-ck  relations, 
neglecting  terms  of  seconi  and  higher  order  in  the  stren/^th  if  the 
shock.  (■■)  and  (c)  are  the  olwlous  extensions  of  the  tvjo  shock  con- 
ditions v;hich  WL.re  used  in  [l]  and  (a)  follows  Mrectl;y^  from  the 
conservation  of  the  tan.rcntial  component  of  m  mentum.   It  may  be 
n-te.]  that  (c)  expresses  the  fact  that  the  shock  vel-ciby  is  the 
mean  :f  the  velocity  of  tho  iN'avelets  on  either  side  of  it;  in  the 
(t,r,z)  space,  it  sh'avs  that  the  shock  surface  bisects  the  angle  be- 
tween the  characteristic  surfaces  meeting  it  fron  either  side. 

It  is  clear  that  the  sh  ck  may  bo  assumed  to  be 

(38)  t  =  a(r,z)  -  7'(r,z)  5^(r,z), 

where  o/a  is  small,  but  a  ;ain  a  crucial  point  Is  that  a'  becomes 

lari:c  compared  to  t-a.   The  component  v.    of  the  particle  velocity 

tan  *' 

is  A,ua^-Awa^  to  the  first  -v\ev,    and  this  must  be  zero  by  (a).   But 

(30)  and  (31)  show  that  the  condition  is  already  satisfiud  by  our 
solution.   The  component  v   is  Ac  u  +  A,,  w  to  first  .-rder,  and  from 

(30),  (31)  this  is  A(a^+a^)XF(7:::1  =XA-^F(r);  hc.nce  (32}  shows  that 

(b)  is  also  satisfied.   The  final  condition  (c)  determines  the 
sh.-ck.   From  (38), 


U  =     •'-    =a!sA+;.^(a  .7"+a  ^  ), 


-r  ^ — z 


since  r/a  is  small;  theref  .-r-. ,  from  (c),  (32)  an  1  tho  abDve  expres- 
sion f'jr  V 


n' 


(39)  a^a'+a^(r  =  lli  X  pf^  . 

r  r  z  z    \\      ,,l\. 

i.  L 

In  addition,  equating  the  values  of  t-a  given  in  (3?)  un  1  {-j>Q)  ,    we 
have 
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(i|0)  a   =  ?F(3:)  -  TT. 

In  principle,  by  eliminating  ~C,  (39)  and  (1+0)  determine  (/ as  a 
function  of  r  and  z.   It  is  more  convenient,  however,  to  find  the 
relation  betvreon  r,  z,  amd  the  value  of— ^  at  the  shock;  this  re- 
lation specif iesl!J'as  a  function  of  r,  z  and  then  (IlO)  gives 
dXv ,z)  .      For  this  purpose,  the  right  hand  side  of  (39)  is  replaced 
from  {}S)    by  ^(a  ^+o.Q)Ft:C),    and  (1+0)  is  substituted  for  cT  to 
give 

(i+1)         |(a^p^+a^p^)F(7D  +  (liP't^-l)(a^^+a^T'.-)  =  0, 
On  multiplyinr:  by  F(7r)  ,  (l+l)  can  be  vjritten 


(1+2)  ^"^va^  ""  °-z^\¥^^^^- l^^'^^i    =  °' 


Since  the  operator  is  the  dfirivaolve  along  the  rays,  the 
term  in  brackets  is  constant  along  each  ray;  but  this  quantity  is 
zero  on  the  initial  wave  front,  therefore 


do) 


|p(r,z)  F^CC)-  j  FfrrMd^C^  =  0, 


at  all  points  on  the  shock.   It  is  of  interest  to  observe  that 
(1+3)  is  exactly  the  formula  which,  vjas  found  for  the  shock  in  the 
problems  involving  only  one  space  variable,  except  in  those  cases, 
of  course,  p  x-jas  a  function  only  of  the  one  variable. 

As  already  rceritioned,  (1+3)  gives  "G  as  a  function  of  r  and  z 
at  the  shock|  from  it,  cTtr,z)  may  be  obtained  by  (!+0)  and  the 
values  of  u,w  and  a  just  behind  the  shock  by  (27),  (2c)  and  (29). 
The  strength  of  the  shock  is  found  sirfiilarly  from 

^^^  P  -  -   P(z)    -    ^2 
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